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Abstract 

Background: Torque Teno Virus (TTV) was the first single stranded circular DNA virus to be discovered that infects 
humans. Although there have been numerous reports regarding the prevalence of TTV from other countries of 
South Asia, there is severe lack of information regarding its prevalence in Pakistan. Thus the present study compiles 
the first indigenous report to comprehensively illustrate the incidence of the virus in uninfected and hepatitis 
infected population from Pakistan. Another aim of the study was to present the sequence of full length TTV 
genome from a local isolate and compare it with the already reported genome sequences from other parts of the 
world. 

Methods: TTV DNA was screened in the serum of 1 16, 100 and 40 HBV infected, HCV infected and uninfected 
individuals respectively. Nearly full length genome of TTV was cloned from a HBV patient. The genome sequence 
was subjected to in-silico analysis using CLC Workbench, ClustalW, ClustalX and TreeView. Statistical analysis was 
carried out in SPSS v17.0. 

Results: Our results report that 89.7%, 90.0% and 92.5% of HBV, HCV patients and healthy control population were 
positive for TTV infection. TTV genome of 3603 bp was also cloned from a local isolate and given the identity of 
TPK01. The TTV genome sequence mentioned in this paper is submitted in the GenBank/EMBL/DDBJ under the 
accession number JN980171. Phylogenetic analysis of TPK01 revealed that the Pakistani isolate has sequence 
similarities with genotype 23 and 22 (Genogroup 2). 

Conclusion: The results of the current study indicate that the high frequency of TTV viremia in Pakistan conforms 
to the reports from other areas of the world, wherever screening of TTV DNA was performed against 5'-UTR of the 
genome. The high sequence diversity among TTV genome sequences and the high frequency of prevalence makes 
it harder to study this virus in cellular systems. 

Keywords: Torque teno virus, TTV prevalence, Co- infection in hepatitis, Pakistani population, TTV genome 



Background 

Torque Teno Virus (TTV) is a circular ssDNA virus [1]. 
It has been classified in a newly characterized genus 
called Anelloviridae. The genome size of TTV varies 
from 3.6-3.8 kb [2]. 

TTV is unique in many ways, it is the first single 
stranded circular DNA virus to be discovered that 
infects humans [3]. Another distinguishing property of 
TTV is its small genome size for a virus that infects 
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humans. Also the nucleotide sequence variation between 
genotypes of TTV is greater than 50% which explains its 
distribution to at least 39 genotypes [4]. Before the dis- 
covery of TTV, this high rate of mutation was never 
observed in a human virus with such a small genome 
size. Moreover TTV genome shows greater variation at 
ORF level than in the UTRs, which is again a distinctive 
trait limited to TTV and TTV-like viruses [5]. 

Reports regarding prevalence of TTV in blood show 
immense disparity, from a viremia rate of 2-23% [6] to 
its frequency in blood rising as high as 90% [7,8]. This 
variation in the rate of positive detection depends on the 
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choice of target site of amplification in PCR, as the 
UTRs of the genome are more conserved as compared 
to the ORFs therefore primers from UTR will cross- 
match a large number of genotypes [9] and hence will 
give a higher rate of detection. 

The conservation of sequence in the UTRs and diver- 
gence in ORFs translates to greater disparity at amino 
acid level than at the nucleotide level. For instance, the 
two TTV isolates SANBAN and TA278 share only 32% 
sequence identity at amino acid level while the nucleo- 
tide sequence of UTRs is 78% identical [10]. Protein ex- 
pression of TTV indicates that three mRNAs are 
encoded from the three overlapping Open Reading 
Frames (ORFs) [11]. The exact function and location of 
these possible protein products is still under study [12]. 

As TTV was discovered in a non-A-G hepatitis patient, 
the prevalence of TTV has been repeatedly examined in 
subjects with hepatocellular carcinoma, Hepatitis B and C 
and at various stages of liver disease. Their results indicate 
a high rate (>90%) of prevalence both in healthy subjects 
and in patients with some form of liver disorder [13-16]. It 
is this omnipresence of TTV and the relative difficulty in 
propagating it in cell cultures, that makes the study of the 
TTV life cycle challenging. 

In relation to hepatitis, some noteworthy observations 
have surfaced. In one study the histological grade score 
was higher in case of HCV-TTV co-infection as com- 
pared to infection with HCV only [17]. There have been 
claims that TTV viral load is a significant factor for 
hepatocellular carcinoma (HCC) in hepatitis C patients, 
irrespective of known risk factors [18]. In another study 
the examination of 25 pediatric cases revealed that chil- 
dren with co-infection of HBV and TTV showed higher 
histological activity index, periportal necrosis and portal 
inflammation scores [19]. 

More than 14 years have passed since the discovery of 
Torque Teno virus but the molecular pathway it follows 
in the host and its pathogenic role, if any, has yet to be 
uncovered. The only report regarding prevalence of the 
virus in Pakistan was made by Prescott et al., in 1999 [20]. 
Thus the present study compiles the first indigenous report 
to comprehensively illustrate the incidence of the virus in 
uninfected and hepatitis infected population from Pakistan. 
The genome of TTV was also cloned and sequenced. The 
sequence was analyzed and annotated through in-silico 
analysis and the phylogenetic relationship with other geno- 
types was determined. 

Results 

Prevalence of TTV detection by PCR 

In this study 116, 100 and 40 samples respectively, of 
HBV, HCV infected and healthy individuals were 
screened with primers from 5' UTR. The results are 
represented in Table 1. Average value of age and ALT 



Table 1 Summary of TTV prevalence in different study 



groups 



Study group 


No. of samples (N) 


TTV positive 
(%age) 


TTV negative 
(%age) 


HBV Patients 


116 


104(89.7%) 


12(10.3%) 


HCV Patients 


100 


90(90.0%) 


10(10.0%) 


Control Population 


40 


37(92.5%) 


3(7.5%) 



levels along with the standard deviations and the male 
to female ratio of the study group are summarized in 
Table 2. 

Statistical Analysis proved that there was no significant 
relationship between factors such as age, gender and 
ALT levels of hepatitis patients and the respective preva- 
lence of TTV in HBV/HCV infected and healthy popula- 
tions. The p-values for each test were above 0.05 except 
in case of Fishers exact test for relationship between 
Age and TTV incidence where p-value was 0.048. Since 
the other two tests (t-test and Pearson s chi-square) gave 
a value >0.05 therefore the result of Fishers exact was 
not sufficient to prove the hypothesis. 

Analysis of cloned TTV genome 

The nearly full length TTV genome sequenced from 
Pakistan was named as TPK01 (GenBank accession 
number, JN980171). TPK01 isolate was 3,603 bp in 
length and had approximately 44 bp missing from the 5' 
UTR and 126 bp missing from the 3' UTR, comprising 
of the GC rich region, and due to primer constraints. 

TPK01 was compared against the nucleotide repository 
of NCBI using Basic Local Alignment Search Tool 
(BLAST) [21] which returned the highest local alignment 
with isolate Kt-08f, genotype 22 (A/C No: AB054647), 90% 
query sequence matched with the subject with a maximum 
nucleotide identity of 77% and E-value of 0.0. 

TPK01 was compared with 106 other full-length 
genomes or full-length cds of TTV in order to obtain align- 
ment for construction of the phylogenetic tree (Figure 1). 

Table 2 Clinical features of the TTV positive and TTV 



negative cases 



Sample group 


Mean 


Gender 


Mean ALT 




age ± SD 


M:F ratio 


Level ± SD 


Hepatitis B 


(N = 1 1 6) 


34.1 ±13.3 


76:40 


90.84 ±47.7 


TTV+ 


34.1 ±12.56 


69:35 


93.35 ±48.7 


TTV- 


34± 19 


7:5 


61.3 ±13.56 


Hepatitis C 


(N = 100) 


38.72 ± 12.2 


51:49 


90.83 ± 40.6 


TTV + 


38.75 zt 1 1.4 


42:48 


90.86 ± 44.2 


TTV- 


38.4 ± 18.5 


7:3 


91.7 ±39.2 
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Phylogenetic analysis was able to conclude that our 
isolate matched most closely to the Genogroup 2 of 
TTV genotypes [22]. Other members of this group in- 
clude, Kt-08f and KMOf from Indonesia [23], ViPi04 
from Italy [24], PMV isolate from United Kingdom 
[25], KAV, TTH18, TTH10, TTH21, TTHllg2, TTH6, 
TTH20 isolates from Germany [26,27]. The NCBI 
Accession numbers of all sequences used in the phylo- 
genetic analysis are indicated in Table 3. Due to 
the high sequence diversity between TTV genome 
sequences reported from around the world, the 
recommended method of classification is to assign 



them in broad genogroups rather than separate 
genotypes [28,29]. 

The sequenced genome was annotated corresponding to 
the highest scoring hit in BLAST results i.e. Genotype 22 
isolate Kt-08f (GenBank Accession No: AB054647) [23]. 
The ORF numbers were assigned according to Genotype 6 
isolate HEL32 (GenBank Accession No: AY666122) [11]. 

The demarcation of ORFs is similar to the Kt-08f isolate 
except for ORF1, which in this case is interrupted by pre- 
mature stop codons at two positions. This could split ORF1 
into three shorter reading frames namely; ORFla, ORFlb, 
ORFlc. The pictorial view of TPK01 is given in Figure 2. 




Figure 1 Phylogenetic analysis of TTV TPK01 based on full-length genome and/or full length cds. The tree was constructed using the NJ 
Algorithm included in the ClustalX (version 2.0.1 1). Bootstrap values of >900 out of 1,000 replicates are shown at the branch points. The tree was 
drawn using Treeview (version 1.6.6). TPK01 (highlighted in yellow) belongs to genogroup 2. The scale bar represents a 10% genetic difference. 
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Table 3 Primers used for screening and cloning of TTV DNA 



S.No Primer name Purpose Sequence Position Expected product size 

1 . T1 F Screening 5-GTA AGT GCA CTT CCG AAT GGC TGA G-3 ; 91-115 1 33 bp 
(NG133) 

2. T1 R Screening 5-AGC CCG AAT TGC CCC TTG AC-3' 204-223 
(NG132) 

3. TTF2 Cloning 5'- TAT GTC GTC CAC TTC CTG GG-3' 50-69 3652 bp 

4. TTR2 Cloning 5'- GAG GAA GGA AGT CGG CCA TTT TG-3' 3684-3707 



Discussion 

Torque Teno Virus (TTV) was first discovered in the 
blood of a post-transfusion non A-G hepatitis patient 
which led to TTV being repeatedly linked as the causal fac- 
tor of cryptogenic hepatitis. There have been numerous 
reports about the TTV viremia in patients of Hepatitis B 
and C from all over the world but there is severe lack of 
data on epidemiology of TTV and its sequence analysis 
from Pakistan. Therefore the present study was designed to 
clarify the rate of prevalence of TTV in hepatitis infected 
population vis-a-vis healthy subjects, outline a relationship 
of TTV prevalence with the incidence of hepatitis and to fi- 
nally analyze the full length genome sequence of TTV from 
Pakistani isolate. 

The findings of this study report a TTV incidence rate of 
89.7%, 90.0%, 92.5% in the Hepatitis B group, Hepatitis C 
group and control population respectively. Our statistical 
analysis refutes the hypothesis that the frequency of TTV 



incidence is influenced by the factors such as age, gender 
and ALT level in hepatitis patients. 

The high prevalence of TTV alludes to the likelihood 
of a latent, possibly lifelong infection. Some of the 
patients in our study were less than 8 years of age. Our 
findings are complemented by other reports of TTV in- 
fection in infants and children [28,30-32] which provide 
evidence that TTV is a mother-to-child transmitted 
virus, as suggested by Yokozaki et al. (1999) rather than 
a blood-borne virus. 

In any epidemiological study of TTV prevalence with 
a PCR based method, the design and target position of 
primers is crucial in determining the frequency of posi- 
tive detection of virus. Due to the high sequence diver- 
sity among TTV genotypes it is improbable that a single 
set of primers will detect all types of TTV. For this rea- 
son the primers used in this study were chosen on the 
basis of extensive review of data pertaining to surveys of 



TATA Signal (38-43) 
5'UTR 



3'UTI 



polyA Signal (2953-2959) 



ETI (nt 302) 
RF2 (302-685) 




MET1I (nt 553) 



,Ary -rich Motif 
Orfla (553-717) 



Splice 01 (nt681) 



ORF3 (2293 2954) 



Splice A2 (nt2492 



Splice A1(nt 2293) 




rflb (802-1251) 



IVR (1362-1610) 



rfic (1651-2757) 



Figure 2 Predicted annotation of TPK01. Splice D, A1, A2, are the splice acceptor and donor sites. The ORF3 will be alternatively spliced to 
code for the putative proteins ORF2/3, ORF2/2, ORF1/2 and 0RF1/1. HVR is the Hypervariable region found in TPK01. The sequence is shown as 
circular for convenience, although the GC-rich tract is missing from the 3-UTR. 
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prevalence of TTV. It has been emphasized by the findings 
of multiple reports that primers from UTRs will have 
greater potential to detect multiple genotypes of TTV ra- 
ther than those designed against the ORFs [33,34]. The rea- 
son being that there is higher sequence conservation in the 
non-coding regions than in the coding regions, where in 
some stretches of the UTRs more than 90% sequence iden- 
tity has been observed [5,35]. 

Based on the above mentioned facts, for the present 
study primer pair NG133/NG132 [33] was selected. A 
thorough analysis of TTV detection strategies by Hu 
et al. suggests that a group of four primer sets have the 
ability to screen all genotypes of TTV and TTV-like 
viruses except genotype 21. Of this recommended array 
of primers, NG133/NG132 is a part [36]. 

The results of the present study are in parallel with 
data published from other parts of the world. For 
example, in the study mentioned above, the positive rate 
of TTV detection was found to be as high as 93% with 
5'-UTR primers but only 10% with primers designed from 
ORF region in the healthy population of Japan. Another 
study which surveyed the TTV infection across different 
geographic locations (Japan, USA, Mayanmar, Korea, 
Nepal, Egypt, Bolivia, Ghana) revealed the frequency of 
positive result between 70-98% from UTR primers, both in 
healthy and disease infected population [16]. The only 
previous report of TTV infection rate from Pakistan was 
made in 1999 in blood donors from Aga Khan University, 
Karachi. The reported prevalence was only 16% with pri- 
mers from the N22 region which belongs to the ORF1 of 
TTV [20]. In separate studies from Russia, Iran, and Hong 
Kong, infection rate in healthy population was 94%, 80% 
and 90% respectively [7,37,38]. In contrast to the above 
reports of high prevalence, TTV was observed in only 2.9% 
of healthy Iranian population emphasizing the fact that 
ORF primers have less conserved sequence which makes 
them incapable of matching multiple genotypes [6]. 

The only studies where the rate of infection in healthy 
population is relatively lower than our findings, when 
similar primers have been used, are from India and UAE 
[39,40]. In these studies nested PCR approach was fol- 
lowed where NG133/NG147 were used in the first round 
while NG134/NG132 amplified in the second round. 
Out of 50 healthy volunteers, only 26 tested positive 
which is a 52% infection rate. The same study shows that 
80.6% of patients with fulminant hepatic failure and 76% 
of patients with Acute Viral Hepatitis tested positive for 
TTV infection. A different research that tested the 
prevalence in patients on hemodialysis and healthy sub- 
jects; the positive infection rate was 83% and 43% re- 
spectively [41]. The study from UAE [40] reported that 
97.9, 95.7% of patients with HBV and HCV tested posi- 
tive for TTV infection which aligns closely with our 
results of 89.7% and 90.0% of HBV and HCV patients 



testing positive for TTV infection. In contrast a mere 
34.9% of healthy population was found to be carrier of 
TTV DNA which differs from our finding of 92.5% of 
control population that was observed as TTV positive. 

The possible reasons for this marked variation in the 
frequency of infection in control population in our study 
and the above cited reviews could be the difference in 
the number of sample population, the fact that the 
prevalent genogroup in these areas is different from 
Pakistan or simply that TTV is more ubiquitous in Pakistan 
than in other countries. 

The TTV isolate cloned in our study named TPK01 has 
sequence similarity spanning over most of the genome with 
the previously reported genotype 22 and 23 isolates Kt-08f 
and Kt-lOf from Indonesia [23]. It is worth mentioning here 
that the previous isolates from Pakistan which were short 
clones of N22 region (200 bp approx.) belonged to geno- 
type la in majority [20]. TPK01 also shows significant simi- 
larity to the PMV isolate from United Kingdom of 
genotype 17 [25]. Stretches of 5'-UTR of the PMV sequence 
are up to 92% identical to TPK01. Other sequences that 
show identity with TPK01 over some stretches of genome 
are ViPi04, isolate from Italy [24], tth6 and tth20 isolates of 
a Hodgkins disease patient from Germany [27]. The entries 
of genome that related most closely to TPK01 were all iso- 
lates belonging to genogroup 2. 

The arrangement of ORFs, promoter and poly-A sites 
is similar to genotype 22 and other members of group 2. 
In TPK01 ntl362-1610 show no similarity with any of 
the TTV genomes assembled in NCBI Nucleotide data- 
base, this region corresponds to the HVR of ORF1 and 
has therefore hypermutated. In a study where TTV iso- 
lates were cloned from Hodgkins patients, the HVR 
extended over ntl425-1693 and showed 80 to 99% variation 
among the isolates that were cloned, whereas the overall 
variation in the sequence of the cloned TTV genomes was 
only 2% [27]. This high rate of variation that translates to 
divergence at the amino acid level could suggest immune- 
evasion by hypermutation [27,42]. Premature stop codons 
have also been observed in ORF1 of our isolate which splits 
the coding sequence into three shorter translation products. 
The occurrence of interrupted ORF1 has been observed in 
other TTV isolates as well. ORF1 splits in two potential 
reading frames in isolates tth7, tt32b8, tth8 while the isolate 
tth22g4 is interrupted twice which translates to three puta- 
tive proteins [27]. TCHN-a, TCHN-cl, TCHN-d2, US32 
also have premature stop codons which can give rise to two 
translation products instead of one [43,44]. 

In one study, subviral molecules could be isolated 
from a Hodgkins lymphoma cell line when it was trans - 
fected with full-length TTV genomes. It was found that 
additional splicing events and rearrangements within the 
genome raised ORFs that could not be explained by the 
parent genomic sequence. These subviral particles varied 
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in sequence length from 172 bp to full length genome. 
Similar kinds of molecules were also isolated from sera 
of pregnant mothers [45]. These findings present the 
question that whether the multiple genotypes and types 
of TTV (Torque Teno Mini Virus, Torque Teno Midi 
Virus) are due to wide intra-genomic rearrangements or 
are actually separately evolved genotypes. 

Conclusion 

The results of this study conform to the previous find- 
ings which have established the high rate of sequence 
variation in TTV genotypes. These variations in se- 
quence are less pronounced in the UTR, but are seen 
extensively in the ORFs. This high level of sequence 
divergence and lack of apparent association with any dis- 
ease poses serious questions about the role of this ever- 
present replicating and mutating form of virus. The 
hardest of all questions is the mode of translation and 
function of TTV proteins. Therefore it is inferred that 
TTV presents itself ubiquitously in human population 
irrespective of any disease association so far. The se- 
quence variation in the genome of TTV can suggest 
rearrangements within the genome giving rise to new 
forms of TTV. 

Methods 

Sample collection 

The study was reviewed and approved by the Ethics 
Committee of Atta-ur-Rahman School of Applied Bios- 
ciences, National University of Sciences and Technology, 
Sector H-12 (NUST), Islamabad, Pakistan. Informed 
consent from all patients and guardians of subjects 
below the age of 18 was obtained before proceeding 
with the study. 

Blood samples were collected from subjects who came 
to Diagnostic Laboratory of Atta-ur-Rahman School of 
Applied Biosciences, NUST, Sector H-12, Islamabad, 
Pakistan for routine testing. 

Confirmation of HBV and HCV infection 

The confirmation for Hepatitis B and C infection was 
done first with serological testing and then the positive 
results were further tested for the presence of HBV 
DNA and HCV RNA. The initial confirmation of infec- 
tion with HBV was carried out with KHB EIA kits for 
anti-HB core, HBe/anti-HBe and HBsAg/anti-HBs 
whereas HCV infection was confirmed with KHB EIA 
kit for the detection of anti-HCV antibodies. 

Further confirmation of HBV and HCV infection was 
carried out with RoboGene HBV DNA Quantification 
Kit which amplified a region of HBV surface gene while 
HCV RNA Quantification Kit was used to test against 
5'-UTR of HCV genome. Standard manufacturer s proto- 
col was followed in all instances. 



Extraction of viral DNA 

DNA was extracted from both HBV/HCV infected and 
uninfected control samples using NucleoSpin Blood 
MACHEREY-NAGEL DNA extraction kit. The DNA 
was eluted in 50 \A of Nuclease-free water. 

Primer designing 
Screening primers 

The primers for screening of TTV DNA (NG133 and 
NG132 as given in Okamoto et al., 1999b, [33]) were 
chosen on basis of their specificity and sensitivity of de- 
tection after extensive review of literature. The primers 
used in the study are listed in Table 4. 

Cloning primers 

Primers were designed manually for cloning of Torque 
Teno virus genome. A total of 23 sequences of full- 
length or nearly full-length TTV genomes retrieved from 
NCBI nucleotide database were aligned in CLC Workbench 
6.1.1 [51], the primers were designed from relatively con- 
served portions. 

Amplification of TTV DNA 

In order to clone the TTV genome, the reaction mixture 
constituted of 5 \A of extracted DNA, 2.5 U each of 
DreamTaq Polymerase (Fermentas EP0702) and Long 
PCR Enzyme Mix (Fermentas K0182), lOpmol each of 
forward and reverse primer, 2 mM ofdNTPs. The PCR 
mixture was preheated at 95°C for 5 minutes, then run 
35 cycles with these parameters: 95°C for 40 seconds, 
60°C for 40 seconds, 68°C for 4 minutes, and the final 
extension was at 68°C for 10 minutes in an AB-Veriti 
96-well Thermocycler. 

For the detection of TTV DNA in HBV/HCV infected 
population and healthy controls the reaction mixture 
constituted of 8 [A of extracted DNA, 1.5 U of Taq Po- 
lymerase, 10 pmol of forward and reverse primer and 
2 mM dNTPs. The PCR mixture was preheated at 95°C 
for 3 minutes, then run 35 cycles with these parameters: 
95°C for 20 seconds, 60°C for 40 seconds, 72°C for 30 
seconds, and the final extension was at 72°C for 7 min- 
utes in a AB-Veriti 96-well Thermocycler. With every 
batch of screening reactions, a positive and negative 
control was used to eliminate the chances of any false 
negative results. A positive control was the serum sam- 
ple which was confirmed for TTV infection after se- 
quencing while the negative control is a serum sample 
repeatedly tested with multiple sets of primers from vari- 
ous parts of the TTV genome but returned negative re- 
sult in all cases (Detail of primers is provided in 
Additional file 1). 

A positive result meant the amplification of a 133 bp 
fragment from T1F/TIR primers. The positive and 
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Table 4 Accession numbers and reference of 106 



sequences used in the phylogenetic tree 



S No. 


NCBI accession 


Isolate 


Reference 




No. 






1 


ARD1 761 D 

ADU I / U I U 


TA278 


niramntn of nl 1 QQfth TR1 
WKdlliUlU cl ill./ I77OU LOJ 


2 


AYR73QRR 

A T OZJ700 


2 h 


NipI pf /i/ 9005 r??i 

INIcI Cl Cff./ ZUUJ LZZJ 


3 


AYR73QRQ 


3 h 


NipI Pf nl 9005 DTI 
INIcI cl Cff./ ZUUj LZZJ 


4 


AY6661 7? 

A I DUO I ZZ 


HEL32 


oin pf ^/ 9005 n 11 

V^IU ct Cll./ ZUUD LI IJ 


J. 


ARD1 761 1 

ADU I / D I J 


Tl KD1 

I UjU I 


Dlramntn pf /?/ 1 QQQr DTI 


6 
O. 


AF761 761 
Arzu I / u I 


PMV 


M^iipt pf ^/ 9000 Dm 

ndiici ct ill./ zuuu l^-~>j 


7 


ARDS4647 


Kt-08f 


Mulinnn Pf /i/ 9001 DTI 

IVIUIJUI1U Ct Cll./ ZUU 1 L' J 


3 


ARDS464R 


Kt-1 Of 


Mulinnn pf 9001 DTI 

IVIUIJUIIU Cl Ml./ ZUU 1 \_/-D\ 


Q 
y. 


AF43^ni4 


KAV 


MpMpi- pf ^/ 9nni r?6i 

nciici ct cii./ zuu i L^-Uj 


1 n 

I u. 


C\ I7Q736D 


BIS8-1 7 


1 Inni iktlickiArJ 
LMl|JUUMollcU 


1 1 


DD1 R7DD6 

L>V^ I O/ UUU 


RM1 CA R 

DIVI I ^_ I O 


1 Inni iktlicho/H 
IsllfJUUMMlcU 


1 2 


ARD64SQS 


CT23F 


p pnn of nl 9009 DR1 
rcliy ct Cll./ zuuz lzoj 


1 3 


ARD64SQ6 

nDUUHJ^U 


CT25F 


Ppnn Pf 9009 DR1 
rcny ct cii./ zuuz l^-OJ 


1 4 


ARD64qQ7 
ADUOtO .7/ 


v_ i jur 


Ppnn pf ^/ 9009 DR1 
rciiy ct cii./ zuuz lzoj 


1 q 
i j. 


ARD646D4 




Ppnn pf ^/ 9009 DR1 
rciiy ct cii./ zuuz lzoj 


1 6 


ARD64SQR 


CT43F 


Ppnn Pf /i/ 9009 DR1 
rcliy ct Cll./ zuuz lzoj 


1 7 


Mmi 4D7S 


CT44F 


Ppnn Pf /i/ 9009 DR1 
rcliy ci cii./ zuuz lzoj 


1 8 


AF1 2291 3 


GH1 


Mu^hahwar pf 1QQQ R1 

IVIUol Idl 1 Wdl Cl Ml./ 1777 \J\ 


1 Q 
I y. 


AF1 77Q1 6 
Ar l zz? l o 


IA1 

JA I 


Fi-kpr pf nl 1 QQQ Dl 
LT Kcl ct Ml./ 1 777 [J] 


zu. 


AF1 77Q1 R 
Ar I ZZ? I o 


IA?h 

J AZU 


Frkpr pf 1 QQQ Dl 
CIK.CI ct Ml./ 1 7^77 L~>J 


21 


AF1 2291 7 


JA4 


Frkpr Pf 1QQQ r 1 

LllVCI Cl Ml./ 1 www L-JJ 


22 


AF1 2291 5 


JA9 


Frkpr Pf 1QQQ r 1 

Lllvci cl Ml./ I www \D\ 


23 


AF1 2291 9 


JA1 0 


Frkpr Pf 1QQQ r 1 

ul l\CI Cl Ml./ I www L-JJ 


24 


AF1 77Q14 3 
Ar I Z-ZSy I t\j 


IA7D 

J AZU 


Frl^pr pf /7/ 1QQQ r 1 
CiKci ct Ml./ 1 www L->J 


ZD. 


ARD64qQQ 


ITD3F 
J I Ujr 


Ppnn Pf nl 9009 T9R1 
rcliy cl Ml./ ZUUZ L^-Oj 


26 


ARD646Dn 

ADUU^-UUU 


JT05F 


Ppnn Pf nl 9009 T7R1 

rcliy cl Ml./ ZUUZ LZOJ 


27 


ARD646D1 

ADUU^UU I 


JT14F 


Ppnn Pf nl 9009 T7R1 

rcliy cl Ml./ ZUUZ L ' J 


ZO. 


ARD646D? 
ADUO^-OUZ 


IT1 QF 
J I I y\ 


Ppnn pf nl 9009 T9R1 
rcliy cl Ml./ ZUUZ L^-Oj 


Ay. 


ARD646D6 


J I jjr 


Ppnn pf nl 9009 T9R1 
rcliy ct Ml./ ZUUZ L^-Oj 


30 


Mr D14D76 


JT34F 


Ppnn Pf nl 9009 T7R1 

rcliy cl Ml./ ZUUZ L ' J 


31 


ARD646D^ 

ADUU^UUj 


JT41 F 


Ppnn Pf nl 9009 T7R1 

rcliy cl Ml./ ZUUZ L ' J 


32 


AYD7646S 

A 1 UZUHUJ 


L01 


1 ill Pf 9009 T4^1 

LIU Cl Ml./ ZUUZ L i-Jj 


DO. 


AYD76466 


I D9 

LUZ 


1 in pf ^/ 9009 T431 

LIU ct Ml./ ZUUZ L^->J 


1A 

DH. 


AF371 ojn 
Arj / I j / u 


I 01 

LUD 


1 iii pf ^/ 9009 T431 

LIU cl Ml./ ZUUZ L^->J 


35 


Mmi4D6Q 


Mrni4D6Q 


VJII|JUlJll3llcU 


36 


AF7QRSRS 
ArzjOJOJ 


P1 C1 


1 Inni iktlicho/H 
UII|JUlJllollcU 


37. 


AR060SQ7 


SAa-01 


Okamntn pf nl 9001 

Valval 1 IV/IU Cl Ml./ £-\J\J 1 L^--'J 


38. 


AB060594 


SAa-10 


Okamoto et o/. f 2001 [23] 


39. 


AB060593 


SAa-38 


Okamoto et al., 2001 [23] 


40. 


AB060592 


SAa-39 


Okamoto et al., 2001 [23] 


41. 


AB060596 


SAf-09 


Okamoto et al., 2001 [23] 


42. 


AB060595 


SAj-30 


Okamoto et al., 2001 [23] 


43. 


AB025946.2 


SAN BAN 


Hijikata et al., 1999 [40] 


44. 


AM712003 


Sle1931 


Leppik et al., 2007 [45] 


45. 


AM712004 


Sle1932 


Leppik et al., 2007 [45] 
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46. 


AM71 1976 


Slel 957 


Leppik et al., 2007 [45] 


47. 


AM712030 


Sle2057 


Leppik et al., 2007 [45] 


48. 


AM712031 


Sle2058 


Leppik et al., 2007 [45] 


49. 


AM712033 


Sle2061 


Leppik et al., 2007 [45] 


50. 


AM712034 


Sle2065 


Leppik et al., 2007 [45] 


51. 


AM712032 


Sle2072 


Leppik et al., 2007 [45] 


52. 


AF247138 


T3PB 


Biagini et al., 2000 


53. 


AF345526 


TCHN-A 


Luo et al., 2002 [44] 


54. 


AF348409 


TCHN-B 


Luo et al., 2002 [44] 


55. 


AF345523 


TCHN-C1 


Luo et al., 2002 [44] 


56. 


AF345527 


TCHN-C2 


Luo et al., 2002 [44] 


57. 


AF345524 


TCHN-D1 


Luo et al., 2002 [44] 


58. 


AF345525 


TCHN-D2 


Luo et al., 2002 [44] 


59. 


AF345522 


TCHN-E 


Luo et al., 2002 [44] 


60. 


AF345528 


TCHN-F 


Luo et al., 2002 [44] 


61. 


AF345521 


TCHN-G1 


Luo et al., 2002 [44] 


62. 


AF345529 


TCHN-G2 


Luo et al., 2002 [44] 


63. 


AB028668 


TJN01 


Ukita et al., 2000 [46] 


64. 


AB028669 


TJN02 


Ukita et al., 2000 [46] 


65. 


AJ620218 


^H3 


Jelcic et al., 2004 [27] 


66. 


AJ620226 


^H4 


Jelcic et al., 2004 [27] 


67. 


AJ620227 


^H5 


Jelcic et al., 2004 [27] 


68. 


AJ620212 


^H6 


Jelcic et al., 2004 [27] 


69. 


AJ620230 


^H7 


Jelcic et al., 2004 [27] 


70. 


AJ620231 


^H8 


Jelcic et al., 2004 [27] 


71. 


AJ620219 


^H9 


Jelcic et al., 2004 [27] 


72. 


AJ620213 


^H10 


Jelcic et al., 2004 [27] 


73. 


AJ620214 


TTH1 1 


Jelcic et al., 2004 [27] 


74. 


AJ620232 


TTH1 3 


Jelcic et al., 2004 [27] 


75. 


AJ620228 


^H14 


Jelcic et al., 2004 [27] 


76. 


AJ620220 


^H16 


Jelcic et a/., 2004 [27] 


77. 


AJ620221 


^H17 


Jelcic et a/./ 2004 [27] 


78. 


AJ620215 


^H18 


Jelcic et a/., 2004 [27] 


79. 


AJ620233 


^H19 


Jelcic et a/./ 2004 [27] 


80. 


AJ620216 


^H20 


Jelcic et a/./ 2004 [27] 


81. 


AJ620217 


^H21 


Jelcic et a/., 2004 [27] 


82. 


AJ620234 


TTH22g4 


Jelcic et a/., 2004 [27] 


83. 


AJ620235 


^H23 


Jelcic et a/./ 2004 [27] 


84. 


AJ620222 


mH25 


Jelcic et a/./ 2004 [27] 


85. 


AJ620223 


mH26 


Jelcic et a/./ 2004 [27] 


86. 


AJ620224 


mH27 


Jelcic et a/./ 2004 [27] 


87. 


AJ620229 


mH29 


Jelcic et a/./ 2004 [27] 


88. 


AJ620225 


TTH31 


Jelcic et a/./ 2004 [27] 


89. 


FR751471 


TTV-HD14i 


de Villiers et al., 201 1 [47] 


90. 


FR751472 


TTV-HD15a 


de Villiers et a/./ 201 1 [47] 


91. 


FR751476 


m-HD16a 


de Villiers et a/./ 201 1 [47] 


92. 


FR751488 


TTV-HD1 7 


de Villiers et a/., 201 1 [47] 
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93. 


FR751489 


^/-HD18a 


de Villiers et al., 201 1 [47] 


94. 


FR751491 


TTV-HD19 


de Villiers et al., 201 1 [47] 


95. 


FR751492 


^/-HD20a 


de Villiers et al., 201 1 [47] 


96. 


FR751498 


TT\/-HD2] 


de Villiers et al., 201 1 [47] 


97. 


FR751499 


TT\/-HD22 


de Villiers et al., 201 1 [47] 


QQ 


rn/ J I jUU 


I I V-nUzod 


ae villiers er ai., zu \ \ |4/j 


99. 


FR751506 


TTV-HD24a 


de Villiers ef al., 201 1 [47] 


100. 


AB038620 


™san-s039 


Takahashi ef al., 2000 [48] 


101. 


AB038621 


™yon-Kc009 


Takahashi et al., 2000 [48] 


102. 


AF247137 


TUPB 


Biagini et al., 2000 [49] 


103. 


AB050448 


TYM9 


Okamoto et al., 2000 [50] 


104. 


AF1 22921 


US32 


Erker et al., 1999 [5] 


105. 


AF1 22920 


US35 


Erker et al., 1999 [5] 


106. 


DQ361268 


ViPi04 


Maggi et al., 2006 [24] 



negative control were also loaded on gel along with the 
gene ruler. 

Cloning of TTV genome 

Due to high sequence diversity within the TTV genome, 
multiple sets of primers were used to amplify the full 
length genome. The TTV genome was successfully amp- 
lified from a HBV positive subject. 

The gel eluted 3.6 kb amplified TTV product was 
cloned in a pCR 2.1 vector using TOPO TA Cloning Kit 
(Catalog no: K4500-01, Invitrogen, USA). The ligated vector 
was then transformed in Top 10 strain of Escherichia colL 
Recombinant plasmids were then prepared from cell cul- 
tures. Initially, the plasmid was confirmed for the presence 
of TTV genome insert with restriction digestion. 

Sequencing 

The cloned TTV genome was sequenced with M13 uni- 
versal primers as well as two sets of internal primers 
(Details provided in Additional file 1) using an ABI 
3730x1 DNA Sequencer (Biobasic Inc Canada). 

In-silico analysis 

Multiple genome sequences available in NCBI database 
were retrieved and aligned with the local Pakistani iso- 
late in CLC Workbench 6.1.1 and ClustalX, version 
2.0.11 [52]. Phylogenetic tree was then constructed using 
Neighbour-Joining Algorithm included in ClustalX. The 
robustness of the tree was tested with bootstrapping 
technique where the tree was replicated 1000 times with 
a seed value of 111. TreeView, version 1.6.6 was used to 
draw the final tree. The protein products from the ORFs 
were predicted in CLC Workbench 6.1.1. Alignment 
with pre-reported proteins was done in CLUSTALW 
and CLC Workbench. Scores for amino acid identity 



and similarity was retrieved from SIAS (Sequences Iden- 
tities and Similarities) server [53]. 

Statistical analysis 

The prevalence data generated with the screening of 
TTV DNA in HBV/HCV and healthy population was 
analyzed with statistical tests carried out in SPSS, ver- 
sion 17.0. The hypothesis to be tested was whether fac- 
tors such as age, gender and ALT levels of hepatitis 
patients affect the rate of prevalence of TTV. Students 
t-test, Pearsons correlation and chi-square/fishers exact 
test were applied wherever relevant. P value <0.05 was 
considered significant. 

Additional file 



Additional file 1: Primers used for establishing negative control 
and sequencing. 
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